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Abstract 


We calculate finite-temperature corrections to the decay rate of a generic neutral (pseudo) scalar 
particle that decays into (pseudo)scalars or fermion-antifermion pairs. The ratio of the finite- 
temperature decay rate to the zero-temperature decay rate is presented. Thermal effects are 
largest in the limit where the decaying particle is nonrelativistic but with a mass well below 
the background temperature, but significant effects are possible even when we relax the former 
assumption. Thermal effects are reduced for the case of nonzero momentum of the decaying 
particle. We discuss cosmological scenarios under which significant finite-temperature corrections 
to the decay rate can be achieved. 
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I. INTRODUCTION 


The cosmological consequences of particles decaying out of thermal equilibrium have 
long been a subject of interest (see, e.g., the early work in Refs. iHsj]). Nearly all studies 
of this kind have neglected the effect of hnite (i.e., nonzero) temperatures on the decay 
rate. This is often a reasonable approximation, depending on the parameters governing 
the decay. However, a few authors have examined thermal effects, with results that are 
scattered throughout the literature. Weldon jsl provided one of the early treatments using 


hnite-temperature field 
appeared in Refs. [10 


heory, the approach we use here. Related calculations subsequently 


15| . Earlier discussions of corrections to the neutron decay rate 


relevant to primordial nucleosynthesis were given in Refs. 


16 


ITII. and corrections to the 


Higgs decay rate into electron-positron pairs can be found in Ref. 


use a different formalism. Later Keil 


18| . but these calculations 


19| and Keil and Kobes 20j reexamined the corrections 


to Higgs decay into e'''e using the real-time formulation of finite-temperature field theory. 
Related calculations for a scalar decaying into fermions were done in Refs, [u, 15|. Recently, 


Gupta and Naya’ 
Czarnecki et al. 


2l| considered corrections to pseudoscalar decay into two photons, and 
22| examined thermal corrections to the decay rate of charged fermions. 


Here we provide a systematic calculation of hnite-temperature corrections for neutral 
decaying particles. Our goal is to provide a more organized and systematic approach to this 
problem in a way that will be useful for researchers in the future. In the next section, we 
provide the formalism for our calculation. In Sec. HI, we examine three cases of interest: (A) 
a (pseudo)scalar decaying into (pseudo)scalars, (B) a pseudoscalar decaying into a fermion- 
antifermion pair, and (C) a scalar decaying into a fermion-antifermion pair. Case (A) and 


case (C) were examined previously in Refs. j^, [l^ and 


19 


3 , 


respectively, but in neither 


case was the enhancement/suppression ratio to the decay rate explicitly studied. Refs. 
[l9, 20| considered a scalar that decays at zero momentum, while our results for case (C) are 
valid for arbitrary momentum for the decaying particle. Case (B) has not been previously 
discussed in the literature. In Sec. IV, we discuss our results and indicate the cosmological 
scenarios to which they are applicable. The most striking effect is the possible enhancement 
of the decay rate for the case of decays into (pseudo)scalars. As we show in Sec. IV, 
an extremely large enhancement is difficult (except for reheating after inflation), but not 
impossible to achieve in the context of the standard cosmological model. We also note that 
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thermal corrections are reduced as the momentum of the decaying particle increases, and 
we provide an explanation for this effect. 


II. DECAY RATES AT FINITE TEMPERATURE 


At zero-temperature, the decay rate 7 /) of a particle with energy Eq can be calculated by 


the Cutkosky rules 


231. This leads to 


7d = 


Im 'Et=o{Eq 


( 2 . 1 ) 


which relates the decay rate to the imaginary part of the self-energy ImSj’=o(T'o) of fhe 
decaying particle and its energy Eq. 

At hnite tem per ature T, the Ciitkosky rules need to be modified. Using the imaginary- 
time formalism 


24 


n 


251], Weldon |9|| showed that for a decaying particle with energy E 


m 


the thermal bath, Eq. (12. Ih is modified into 


Tn±Tj = 


ImS(E) 


( 2 . 2 ) 


where T^ is the finite-temperature decay rate, “-I-” and correspond to a decaying fermion 
and boson respectively, and T/ is the inverse decay rate of the particles resulting from the 


decay. Up to one-loop calculation, t 
who used the real-time formulation 


lis result was confirmed by Kobes and Semenoff 


24 


26 ] 


25(1 . If the unstable particle decays in a thermal 


bath that is abundant in its decay products, the decay products would have the proba¬ 
bility to recombine in the thermal bath, and ImS(E') accounts for both of the decay and 
recombination processes. 

Weldon 9(| also showed that regardless of whether the decaying particle is a fermion or 
boson, the ratio of T^ to T/ is a universal function of U, namely 

Td 


= exp (PE), 

^ I 

with (3 = 1/T. This allows us to derive the decay rate at finite temperature 

1 




ImS(U) \ 


(2.3) 


(2.4) 


1 ± e-is® V E 

where again “-I-” and correspond to a decaying fermion and boson respectively. 

In this paper, we are interested in an unstable particle that is out of equilibrium. We 
assume that the finite-temperature corrections to the mass of the decaying particle are 
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FIG. 1: (Left) Self-energy for <I> with a (pseudo)scalar (p loop. (Right) Self-energy for $ with a 
fermion / loop. 


negligible compared to its mass in the vacuum. So we can approximate E as Eq. As we shall 
see, the imaginary part of the self-energy of the decaying particle can generally be written 
as a linear combination of zero-temperature and hnite-temperature contributions, and with 
the approximation E Eq^ we can write ImS(i?) ImS 7 ’=o(-E'o) + IniS 7 ’^o(-E'o)- We can 
then define the ratio 


Id 


(2,5) 


which characterizes the missing factor we would encounter if we blindly use the zero- 
temperature decay rate yz) in a thermal bath. The calculation of R (generalized to arbitrary 
momentum for the decaying particle) for the cases of interest is the main goal of this paper 


(and the results that extend this work beyond that of |9|, ll5|, ll9|, l20| ) 


III. SPECIFIC PARTICLE DECAY RATES AT FINITE TEMPERATURE 

We illustrate our study by considering three simple models: (A) a (pseudo)scalar decaying 
into (pseudo)scalars, (B) a pseudoscalar decaying into a fermion-antifermion pair, and (C) 
a scalar decaying into a fermion-antifermion pair. In particular, we study the ratio R = ^ 
and investigate how it changes with temperature. All of the calculations are done under 
the imaginary-time formalism. This formalism has the advantage that perturbation theory 
can still be organized into a diagrammatic expansion with the same vertices as at zero 
temperature. 
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A. (Pseudo)scalar Decaying into (pseudo)scalars 


We consider the model in which a (psendo)scalar $ can decay into a pair of identical 
(pseudo)scalars (jxj). The interaction operator responsible for this process is: 


C-mt = 9 4 ' 


(3,1) 


This model is relevant to several cases of interest. For instance, <F conld be the Standard 
Model (SM) Higgs decaying into a pair of scalar dark matter particles |27|, or conversely, a 
scalar dark matter particle decaying into a pair of SM Higgs. Alternately, $ could be the SM 
Higgs decaying into a pair of light CP-odd scalars in the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) (28|, or the heavier CP-odd Higgs decaying into two lighter CP- 
even Higgs in the two-Higgs-doublet-model (2HDM) Finally, <h could be the SM Higgs 


decaying into a pair of pseudo-goldstone bosons, which has been proposed by Weinberg 
to explain the fractional effective number of neutrinos hinted by Ref. 31|. 


^1 


To apply the Cutkosky rules, we need to calculate the self-energy of <F as shown in Fig. 
[T] (left) and then put the (j) particles on their mass-shell. Based on the calculations (of the 
imaginary part of the self-energy) in Appendix Al, Eq. (12. 4 h and Eq. (12.Ih . we obtain the 
rates for the decay <F ^ 0 0 at both zero and hnite temperatures: 


7d = 


D 



4 

-Jp 0[M|-4m2] 


1 - 


Id + 


g^T 
4:71 Eok 


In 


1 - 
1 — 


0[M| -4ml 


(3.2) 

(3.3) 


R. 




1 - 


1 + 


2T 


In 


7 / H 4 m? 


1 - 

1 — 1'^ 


0[M| 4m|], (3.4) 


where 0 is the Heaviside step function and 


En ± k \ 1 


Eo = JP + M|, u;± = 


4mj, 

_0 

M|" 


(3.6) 


with k being the momentum of the <F particle. 

Typically, would receive hnite-temperature corrections which go like ^ T where ^ 
is a perturbatively small constant. For ^ < 0.01 and T < 50M,j,, the hnite-temperature 
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FIG. 2: The plot of against T/M^ for a nonrelativistic assuming 4m^ ^ M|. The solid 

(red), dashed (blue), dotted (purple) and dot-dashed (green) lines correspond to the parameters 
k/M^ = 0.001, 0.25, 0.5, 1 respectively. 

corrections to are negligible compared to and therefore can be ignored in the quantity 

ml/Ml 

For the calculations in Appendix Al, we have used the thermal propagators for the 0 
particles and so they are required to be in thermal equilibrium with the thermal bath. It 
is precisely these thermalized (j) particles that induce hnite-temperature corrections to the 
decay rate of $ ^ 0 0. To ensure that there is a signihcant abundance of the (j) particles in 
the thermal bath, we also require that T > mff). 

Now consider an out-of-equilibrium (pseudo)scalar <F which decays into a pair of identical 
(pseudo)scalars (jxp. We plot as a function of T/M^ in Fig. [2l Fig. [3] and Fig. 01 

taking <C M|. This last assumption is not essential and is only used to simplify the 
plots; we have verihed that relaxing this assumption gives similar results. 

It is clear that for a nonrelativistic <F with k/Mq, < 1, the enhancement factor can 

be as large as 10^ (see Fig. ED. For a slightly relativistic <F, the enhancement factor can be 
at least 10^ for T > 20M$ and can reach 10^ for T ~ 50M$ (see Fig. [3D. Even for a highly 
relativistic <h, the enhancement factor can still reach 10 or higher for T ~ 50M,j, (see Fig. 
ED- We thus conclude that the condition T/Mq> ^ 1 is the key for large thermal effects on 
the decay rate. The magnitude of the thermal effects increases signihcantly when we move 
from the relativistic limit to the nonrelativistic limit. 

In the discussions above, we have considered T as large as 50M$. We have not taken into 
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FIG. 3: The plot of against T/M^ for a slightly relativistic ‘h, assuming 4m^ <C M|. The 

solid (red) and dashed (blue) lines correspond to the parameters k/M<^ = 5, 10 respectively. 
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FIG. 4: The plot of against T/M^ for a highly relativistic assuming 4m0 ^ M|. 

The dotted (purple) and dot-dashed (green) lines correspond to the parameters k/M^ = 50, 100 
respectively. 

account the possible finite-temperature correction to M|. With a more accurate calculation, 
M| at finite temperature T should take the form 

Ml{T) = Ml + AMliT), (3.6) 

where AMl{T) arises from the real part of the self-energy for <|). The quantity AMl{T) 



0 10 20 30 40 50 


r/Mo 
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represents the finite-temperature correction to M|, and it is computed in Appendix A2: 


.1 


^2 _ u 


24 




2 j^2 

2 


m 


(3.7) 


0 


Note that the real part of the self-energy for $, from which AM^{T) is extracted, was not 


computed in Ref. |9| but discussed in Refs. Idl. 1,32 1. 


For a nonrelativistic $ with k/Mq, < 1, we have ~ ~ M|, and so AM|(T) 

For consistency, we require AM|(T) < M|. In order to obtain 2> 1 in Fig. 

[21 we have taken Mq, < T/5, which gives: 


9lll< 

24 ~ 




Isl $ . 


(3.8) 


Therefore, is the consistency condition that allows one to neglect the finite- 

temperature correction to M| for the range of T/M^ values shown in Fig. |2l 

In contrast, for a relativistic <F with k/Mq> > 1, we have u'^ — k'^ ~ 0, and so AM|(T) 0. 
Thus, one can always neglect the finite-temperature correction to M| for the range of T/M^ 
values shown in Fig. [2] and Fig. jH There is no analogous upper bound on \g\. 

A final remark is in order. When the temperature is sufficiently higher than the masses of 
the particles under consideration, naive perturbation theory may break down, especially for 
soft external momenta (fc -C T). One could correct this by using a resummed perturbative 
expansion 33|. See Ref. for recent work that discusses this issue. 


B. Pseudoscalar Decaying into a Fermion-Antifermion Pair 


We consider the model in which a pseudoscalar $ can decay into a fermion-antifermion 
pair //. The interaction operator responsible for this process is: 


7:int ^ *A/7^'5/, 


(3.9) 


This model could be relevant in several cases. For instance, <F could be a pseudoscalar dark 
matter candidate, which is the neutral component of an SU(2) multiplet, decaying into SM 


fermions 


34| . Alternately, <h could be the Majoron decaying into Majorana neutrinos 35| . 


^ Notice that in the limit mq, —>■ 0, the quantity in Eq. (IA30|) becomes divergent. This is a manifestation 
of the infrared divergence due to massless particles at finite temperature. The approximation made in 
Eq. (IA27I) may no longer be consistent. In this case, a simple analytical form for AM|(T) may not be 


available. 


















To apply the Cutkosky rules, we need to calculate the self-energy of as shown in Fig. 
[1] (right) and then put the fermions / and / on their mass-shell. Based on the calculations 
in Appendix B, Eq. fl2.4p and Eq. fl2.ip . we obtain the rates for the decay $—)>// at both 
zero and finite temperatures: 


Id = 

Fz, ^ 


A^Ml / 4m2 

SttEoV Ml 


1 

1 - 


Id + 


0[M|- 

X^MlT 
4 : 7 t Eok 





In 


1 + \ 
l + e-i?-/T j 


0[M|-4mJ] , 


(3.10) 

(3.11) 


^ ^ 1 

~ ^ _ q-EoIt 


2T 



In 


1 + A 
1 + e-®-/'^ j 


0[M|-4mJ], (3.12) 


where 




Eq i k 



2 


(3.13) 


Similarly, would receive finite-temperature corrections which go like T where 
is a perturbatively small constant. For < 0.01 and T < 50M$, the hnite-temperature 
corrections to m/ are negligible compared to and therefore can be ignored in the quantity 
mj/Ml. 

For the calculations in Appendix B, we have used the thermal propagators for the fermions 
/ and / and so they are required to be in thermal equilibrium with the thermal bath. It is 
precisely these thermalized fermions / and / that induce finite-temperature corrections to 
the decay rate of $ —)■ / /. To ensure that there is a significant abundance of the fermions 
/ and / in the thermal bath, we also require that T > mj. 

Now consider an out-of-equilibrium pseudoscalar $ that decays into a fermion-antifermion 
pair, f f. We plot against T/Mq, in Fig. [5l As we can see, the suppression factor 

Rq^ff does not vary much when T/Mq increases from 1 to 50. When one moves from 
the nonrelativistic limit to the relativistic limit, the thermal effects decrease and hence the 
suppression factor increases (less suppressed). 
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FIG. 5: The plot of against T/M^^ assuming Am'j <C M|. The solid (red), dashed 

(blue), dotted (purple) and dot-dashed (green) lines correspond to the parameters k/M^ = 
0.001, 1, 10, 100 respectively. 

C. Scalar Decaying into a Fermion-Antifermion Pair 


We consider the model in which a scalar $ can decay into a fermion-antifermion pair, 
//. The interaction operator responsible for this process is: 




(3,14) 


For instance, $ could be a scalar dark matter candidate, which is the neutral component 
of an SU(2) multiplet, decaying into SM fermions [3^. Besides, $ could be the SM Higgs 
decaying into SM fermions. 

As in the pseudoscalar case, in order to apply the Cutkosky rules, we need to calculate 
the self-energy of <F as shown in Fig. [T](right) and then put the fermions / and / on their 
mass-shell. Based on the calculations in Appendix C, Eq. (12. 4 h and Eq. (I2.1jl . we obtain 
the rates for the decay <h ^ / / at both zero and hnite temperatures: 

3/2 

0[M|-4m^], (3.15) 


= 


y‘^M\ 


F 


_$ 

1 ttEq 

1 


ArrA 


D 


1 - 


Id + 


y^MlT 
A'k E^k 


1 - 


4 m 


Ml 


In 


. _ jrp 

ryyr I 0[M| - 4mj] 


1 + e 


-E-/ 


(3.16) 
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1 




1 + 
1 + 


0[M|-4mJ], (3.17) 


< 1 -^// 1 _ ^-EqIT 


1 + 


k 



where is given by fl3.13p . Again, ruf would receive finite-temperature corrections which 
go like ^"T where is a perturbatively small constant. For < 0.01 and T < 50M$, the 
finite-temperature corrections to ruf are negligible compared to and therefore can be 
ignored in the quantity Moreover, we require that /, / are thermalized and T > ruf 

for reasons explained in the pseudoscalar case. 

We find that is identical to R^^/f and so we can just refer to Fig. [5] for the 

behavior of against T/Mi^. The conclusion is similar to the pseudoscalar case. 

IV. DISCUSSION AND COSMOLOGICAL IMPLICATIONS 

The results presented here are in broad agreement with our intuition from statistical 
mechanics. For decays into fermions (III.B. and III.C.), the result of finite-temperature 
effects is a suppression of the decay rate for T/M$ S> 1 resulting from Pauli blocking. 
Conversely, for decays into (pseudo)scalars (UFA.), one sees significant enhancement from 
stimulated decays when T/M$ 3> 1. 

Our results show that thermal corrections are reduced for nonzero momentum of the 
decaying particle; in all of the scenarios we explored, the finite-temperature effects decrease 
as k/Mq, increases. Note that this is not due to Lorentz suppression of the decay rate, as the 
ratio R defined in Eq. fl2.5p includes the same Lorentz factor in both the numerator and the 
denominator. To understand this effect, suppose that /c/M$ is large, and transform to the 
rest frame of the decaying particle. In this frame, the thermal background has a large net 
nonzero mean momentum. But the particles in the thermal bath that interact to produce 
a reverse decay must have zero total momentum, a result that becomes more difficult to 
achieve as the thermal background is boosted to higher and higher momentum. 

When are these results relevant for cosmology? The only cosmologically-relevant decay 
process known to occur with certainty is the decay of free neutrons into protons during 
primordial nucleosynthesis, which occurs at a temperature T ~ 0.1 MeV. In this case. 


■ tron ^ 1, so we expect thermal corrections to be very small, as they indeed are 
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Now consider more hypothetical scenarios. As we have seen, a large change in the decay 
rate occurs only for T/M$ ^ 1. For a particle with a standard thermal history that drops 
out of equilibrium when it is nonrelativistic, we automatically have T/M$ <C 1, so if this 
particle subsequently decays, thermal corrections to the decay rate will be negligible (for 
this and other scenarios discussed here, see, e.g. Ref. |36|h 


On the other hand, if the particle drops out of equilibrium while still relativistic, we 
would have T/M$ ^ 1 when this decoupling occurs. However, in this case /c ~ T at all 
later times, so that ~ T/M^ S> 1. Thus, if the particle decays when TjMq, ::§> 1, it is 

still relativistic at decay {k/Mq, :$> 1). In this scenario, decay into (pseudo)scalars can still 
produce an enhancement of 0(10) (see Fig. 4), but not the O(IO^) enhancement in Fig. 2. 
To achieve the latter requires the transfer of entropy into the thermal background so that 
T/M$ 1 when /c/M$ 1. Some entropy transfer occurs in the standard cosmological 

model when particles that are in thermal equilibrium become nonrelativistic 3^. A larger 
effect can occur in nonstandard scenarios when nonrelativistic particles come to dominate 
the energy density of the thermal background and then decay out of equilibrium In 
either case, the thermal background will be heated so that k < T, and one could then have 
a decaying particle with T/M^ 3> 1 and k/Mq, 1. (This loophole is in principle possible 
even when <F decouples while nonrelativistic, but it would require an enormous entropy 
release in this case). 

A third possibility is a nonthermal production mechanism for the decaying particle in 
question. For example, axions (or axion-like particles) produced by the misalignment mech¬ 
anism are “born” with T/Ma ^ 1 and k/Ma <C 1. 

One possible cosmological scenario for which thermal effects might be signihcant is re¬ 
heating after inflation. Once the inflatons start to decay, the decaying products may form 
a dense plasma which back-reacts on the inflaton decay It is possible that T/M$ ::§> 1 
is satished during this period and our results apply. The effect of this dense plasma on 
the thermal history of the universe was investigated in Ref. 38|. Our results may also 
be relevant to the fate of flat directions after reheating. It was pointed out that thermal 
corrections could be signihcant in this context [h, 39|. 

Thus, while the conditions necessary for thermal corrections to produce an extremely 
large change in the decay rate are somewhat unusual (except for reheating after inhation), 
they are not impossible to achieve in the context of our current cosmological model. Of 
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course, our results are also valid in the case of smaller corrections to the decay rate, which 
are easier to achieve. 
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Appendix A: <1) 


1. Decay Rate: Imaginary Part of the Self-Energy 


The treatment in this subsection is similar to that of 
the held $ in the Matsubara representation is given by 

1 


10| . The one-loop self-energy of 


S( l^ni ^ ) 2 ^ J ^27j-^3 ^ ^ P ) ^(pi, T ^ni P T ^ ) 


(Al) 


where Um = 27rm//5 and = 27r?7,//3, with m, n = 0, ±1, ±2, ..., are the bosonic 
Matsubara frequencies. The Matsubara propagators are written in the following dispersive 
form: 

Pi(Po, P) 


G^{uJrn,P) = J dpo 
G(l){uJrn + P + k) = J dqQ 

where the spectral densities are 


PO * 

p2{qo,P + k) 
Qo ^ ^ 


PiiPo,p) = TT—[<5(^0-^^!)-^(Po + ^^i)], = Ji^ + ml , 

Z CJi V 


:uji 

1 


(A2) 

(A3) 

(A4) 


P 2 iqo,P + k) = - — [Siqo-uj 2 )-Siqo + U 2 )], UJ 2 = J{p + k^ + ml . (A5) 

ZUl2 * 

This representation allows us to carry out the sum over the Matsubara frequencies ojm in a 
rather straightforward manner 2J, [2^ : 




^ PO-lUm qo-lUJm-ll2n 

<^m 


nsiPo) - nBjqo) 
qo-Po-ii^n 


(A6) 
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where nsioj) = is the Bose-Einstein distribution function. The resulting self-energy 

can now be written in the dispersive form: 


E(z/„, k) 




ImE( cu, k ) 
u — ifn 


(A7) 


with ImS(ci;, k) being the imaginary part of the self-energy 


ImS(a;, k ) 

= -2 7r/ 


d?p 

(27r)3 


dpodqQ [nB^Po) - ^^(go) ] Pi(Po, p) p2iqo: P + k) 6{uj - qo + po) , 

(A8) 


The retarded self-energy is dehned by the analytic continuation: 


Sret( ko, k) = = -i ko - e, k) = 


du 


ImSret(w, k) 


71 


(A9) 


u — ko + ie 

Integrating over dpo and dqo, using the identity hb^—oj) = —( l+UBioj)) and performing the 
transformationp —p—k in all the integrals involving nB{oJ 2 ), we can write ImSret( OJ, k) = 
(To -|- (Tt where 


(To = 

(T'p = 


9^ . r \ f d^P 


16 TT^ 
-.2 


sign(a;) 


OJl OJ 2 


5{ |a;| - (Ui - 072 ), 


r d^p 

sign(a;) / - UBioJi) S{\uj\- 001 - 002 )- 


8 tt' 


OOi 002 


(AlO) 

(All) 


Obviously, (Tq represents the zero-temperature contribution while (T^ gives the finite- 
temperature correction. Notice that there were some possible terms involving 5{oo+ooi—oo2 ) 
and 5{oo — 001 + 002 ) in ImSret( Ci;, k ), but they are kinematically forbidden. To proceed, let 
O = cji and z = 002 - Then, we have 

^2 poo nz'^ 


O'O + — 


9 


Sti k 


sign(c(;) / [ 1-|-2 77 , 5 ( 0 ) ] (iO / 6{\oo\ — Q — z) dz , 


(A12) 




where z^ are given by 


z"^ = \J {p^ k)"^ -|- = ^jo? ±2k d Q? — + kP. 


(A13) 


For the integral to be non-vanishing, we require that 


.2 < Z = \u\ — kl < Z'^ . 


(A14) 
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Squaring both sides twice properly, these two inequalities can be cast into the condition 
f{Q) < 0 where 

/(hi) = A{\oj\^ — — A\oj\{ \oj\^ — k‘^ ) kl + {\uj\‘^ — k‘^Y + 4 . (A15) 


Notice that the graph of /(f^) against hi represents a conic with a positive y-intercept. 
Solving f{kl) = 0 for hi, we obtain two solutions: 


\u 


= 


(Icup — /c^ ) ± /c (IcuP — k'^Y — A{ IcuP — k‘^)m?^ 






(A16) 


There are two possibilities: (i) |a;p — /c^ > 0, (ii) k"^ — \oj\^ > 0. For /c^ — |a;p > 0, the 
graph with f{kl) against hi shows that the condition (IA14D can be satished only if hi > u:~ 
but algebraic calculations indicate that |a;| — u:~ < 0. Hence, the condition flA14p cannot 
satisfied and this solution should be discarded. 

For Icnp — > 0, a detailed analysis of /(H), and |ci;| — H as functions of H reveals 

that that condition flA14ll can always be satisfied as far as u:~ < H < and |a;| > 
+ m^. For the discriminant in oj^ to be positive, we require |a;| > 
or |a;| < k. Since > k, we can only choose |a;| > y^A:^ + 4 


0 ' 


As a result, using the integration formula J In (1 — e we conclude that 

ImSi.et(n;, k) = ao + <Jt with 


do = 


CT-r — 


{uj~^—uj ) sign(a;) 0[ |a;|^ — A:^ — 4m^], 


Sti k 

,2 / i_ 
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All k (3 


In 


1 — e~^‘^ 


sign(a;) Q[\u\‘^ — k'^ — 4m|], 


(A17) 

(A18) 


where can now be safely simplified to become 


= 


Icnl ± A: \/ 1 — 


4m^ 


(A19) 


2. Dispersion Relation: Real Part of the Self-Energy 


The real part of the self-energy is given by 


ReS(i/„, k) = 2g^ 


d^p 1 




(27r)3 p ^ p2 + ml + ul (^+ k)2 + rnl + {um + i^n)^ 


. (A20) 
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To proceed, we introduce the Schwinger parameters 




-ai ) 


dai e 


dao e““^ [ (P+^)^+"l| + ( ‘-^m + Un )2 ] 


(A 21 ) 

(A 22 ) 


{p +k)'^+ ml+{um +i^n)^ Jo 

After completing squares, the d^ p integrals become Gaussian and can be easily evaluated 


to give 


ReS(z/„,fc) = 


8 TT^ 


I doii I d(y.2 - 2 

'0 Jo ( tti + Q!2 ) 




= al+a2 


na2 i 1 


cxi T cr 2 4 TT tti T tt 2 
where k'^ = k"^ + The Jacobi theta function 'd^z, r) is dehned as 


^) = X] 


2 2 TT m z-\-i TT m? r 


, (A23) 


(A24) 


^ • 2 / 

and we have used the identity -d^z, r) = (—/•>■ 'd{z/T, — 1 /r). 
Let a = «! + 02 and x = (ai — a 2 )/a. Then, we obtain 

9 /•CXD J 

ReS(z/„,fc) = 


16 TT^ ./n a 




dx o(l -a:), 


'-1 


i/3^ 
4:71 a 


To perform the integration over dx, we can use the formula 


/ dx --— 

1-1 2y/B 


2 -t\/7r ^-Ay4B 


(A25) 


.R + A/2\ r.B-A/2\ 

erj i -—— + erf ^-;=— ) 


, (A26) 


/By \ Vb J 

where erf{z) = fj dt is the error function. In this problem, A = innm and 
B = ^ak"^. The leading contribution of the integral flA25ll comes from the region a ~ 0. 
Near this region, we have 

i k% 


.B + A/2\ , ^f.B-A/2 

erf i -=— + erf i 


Vb J \ Vb 

Meanwhile, ReS(z/„, k) can be written as 


E „ T o fci ^ ak'i. 


2 77,2 ^2 


g4“"'iS g 


.(A27) 


+ 00 


+ 00 


ReS(z/„,fc)= ^ J^ = /q+ ^ J, 


m=^0 


(A28) 
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The quantity Jq corresponds to the zero-temperature contribution and we assume that it 
has already been combined with the bare mass-squared of <h to give M|. Therefore, the 
mass-squared of $ at hnite-temperature T takes the form M^{T) = M| -|- AM|(T) with 
AM|(T) being the finite-temperature corrections: 

+ CXD 

AMl{T)= (A29) 


Upon some simplihcations, we get 


‘■rn^O 


9^ kl 
32 TT^ 77,2 ^2 


da e 


2 1 8 ^ 
—am^ — — — 

<p a \ 4 


(A30) 


which is an even function of m and so AM^{T) = 2 X]m=i ^rn- We can perform the remaining 
integration using the identity 


da e 


-aC-^D 


= 2 


Ki{2VCD), 


(A31) 


where Ki{z) is the modihed Bessel function of second kind. 

Applying the analytic continuation: —ioj — e, we hnd 


H-CX) 


AM|(T) 


„2 , ,2 1.2 

9 CU — ft ^ / o 

ir ~2 - —Ki{m(3 

^ m=l 


niA 


(A32) 


Since Ki{z) ~ for 2 ; 1, it is obvious that ^Ki{m(3m^) will be exponentially 

suppressed if ml3m^ ^ 1. On the other hand, Ki{z) ~ ^ for 2 ; -C 1. Thus, the dominant 
contribution of — Ki{m BmA goes like — —i—. Using ^ we obtain 

jyi L\ r-' cpj Q Yn mpm^ ® Z—/m=l 6 ’ 


AM|(T) ^ 


24 


a;2 - e T2 

a ;2 


(A33) 


Appendix B: Pseudoscalar 4> —)• // 


The one-loop self-energy of the field <h in the Matsubara representation is given by 

d^p 1 


/c) = -A^ 


(27r)^ /3 


^ Tr Gf{um, p)'J^ Gf{um +^n, P+ k)'y^ 


(Bl) 


where = 2 7r (m -t- ^)//3 and = 2 vr (tt, -|- i)//3, with 771,77 = 0, ±1, ±2, ..., are the 
fermionic Matsubara frequencies. It is convenient to write the Matsubara propagators in 


17 













the dispersive form: 


Gf{Um,p) = J dpo 
Gf{Um + iyn,P + k)= / dqo 


pijpo, p) 

PO i 

P2(go, p + k) 
QO ^ 


( TPo-l-P + ^fr.! N e/ , M 

Pi(Po,p) = -^-^[(5 (po-^^i)-<5(po + ^^i)J, 

Z LUi 


(B2) 

(B3) 

(B4) 


P2(?o, P + *:) = 


7° (^0 - 7 ■ (P + fc) + 

2uJ2 


[S{qo-U2)-S{qQ + uj2)], (B5) 

(B6) 


= ^Jp^ + rri^f, a;2 = Y (p + /c)^ + 

This representation allows us to carry out the sum over the Matsubara frequencies ujm in a 
rather straightforward manner 2J, l25(| : 

1 -sp 1 1 _ npipo) - npiqo) 

fd 


(B7) 


PO-tUJm qO-^(^m-^JZn qo - PO - ^ 

'^m 

where npioj) = is the Fermi-Dirac distribution function. The self-energy can be 

written in the dispersive form: 


E( 1/^, k) = -- 


dco 


71 


ImS( cj, k 
u — i Ur, 


(B8) 


where ImS(a;, k) is the imaginary part of the self-energy given by 
ImE(a;, k ) 


= ttA^ 


d?p 


dpo<*90 IKf(Po) - ’>f(9o)l Ti' ( Pi(po,p)7''P2(9o, P + ^7^ <*(i*) - 9o + Po), 

(B9) 


We can then proceed by using Tr(l) = 4 and ^(7^7*^) = giving 


Tr 


(7°Po - 7- PT nr) 7 ^ ( 7°- 7 ' (P + ^) + nr ) 7 


= -4 ( po <?o - P ■ (P + ^) - nr^ 


(BIO) 


The retarded self-energy is defined by the same analytic continuation as in Eq. 


Similarly, integrating over dpo and dqo, using the identity np{—u) = 1 — npioj) and per¬ 
forming the transformation p —p — k in all the integrals involving 7ri?(a;2), we can write 
ImEi.et(CiT, k) = ao + (Tt where 




(To = 

CTp = 


8 vr^ 

A2 


sign (or) 


d^p 


dTT^ 


sign(a;) 


OJl 0 J 2 
d?p 


OJ 1 OJ 2 + p ■ {p + k) + mj \ S( Icnl - cji -U 2 ), (Bll) 


CJl CJ2 


npicui) f cJi 0*2 + P • {p + k) +mj j 6{ \u\ - cui - UJ 2 
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Again, ao represents the zero-temperature contribution while ar gives the finite- 
temperature correction. Also, there were some possible terms in ImSret(i^, k) involving 
5{uj + UJ 1 —UJ 2 ) and 5{uj — UJ 1 + UJ 2 ) which are kinematically forbidden. To proceed, we again 
let n = cji and 2 = 002 - Then, we have 5{\uj\—uji — uj 2 ) = (5( |a;| — hi — 2;), and we can make 
the following simplification: 

ujiUJ 2 + p-{p + k)+m) = -^^-, (B13) 


using the constraint (5( |a;| — hi — 2;). This leads to 

^2 /-oo 


<^0 + — 


8 TT /C 


(Icup —) sign(a;) / [l — 2nF{^l)]d^l / 6{\uj\—^l — z)dz, (B14) 


where 2;^ are given by Eq. flAlSp with replaced by m/. We can then follow the similar 
kinematical arguments in Appendix A to facilitate the integrations over both of dQ and dz. 

As a result, using the integration formula f In (1 -1- we conclude that 

ImEi.et(n;, k) = ao + ar with 
A2 


(To = 

(Tj' = 


8 TT fc 


(Icup — — E ) sign(a;) 0[ [cup — /c^ — 4mj], 

1 + e-^®^ 


(B15) 


Auk jd 

where E^ is given by 


(leap -k"^) In 




sign(c(;) 0[ jeuf —/c^ — 4my ], (B16) 


E^ = 


|ci;| ± k \ 1 




(B17) 


Appendix C: Scalar $—)•// 

Similar to Appendix B, the one-loop self-energy of the field $ in the Matsubara repre¬ 
sentation is given by 

where Um = 2 vr (m A- \)/fd and = 2 7r {n + \)/P, with m , n = 0, ±1, ±2, ..., are the 
fermionic Matsubara frequencies. 


^ f ( ^mi P ) ^ f ( T ^ni P d~ k') 


(Cl) 
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Following the similar steps and tricks as in Appendix B, we obtain ImSi.et( k ) = ao+aT 
with 

(To =- - ( \u}\^ — — 4mj ) ( — E ) sign(c(;) 0[ \bj\^ — /c^ — 4m^] , (C2) 

8 TT ru 

0-2 / 1 + \ 

^ (|a;|2-/c2-4mj) In ^ j sign(a;) 0[ |a;|2 -/c^ - 4mJ] (C3) 

where E^ is given by (1B17F 
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